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Abstract. We investigate the possibility of identifying 
massive objects (lenses) in the Galactic Center region 
(GC) by means of pulsar timing. The well known intensity 
change due to microlensing is found to be less important. 
For typical stellar masses, the flux magnification can be 
significant only if the lens passes very close to the pulsar- 
Earth axis. We show that in the case of a pulsar the time 
varying travel-time delay, which is observable because of 
On the pulsating nature of pulsar radiation, is a much more 
powerful tool to investigate mass distributions in the GC. 
' We find that a travel-time delay is measurable even for 
rather large distances between the lens and the pulsar- 
^ Earth axis. 

^— ( The time varying travel-time delay can be used to de- 
termine the mass of the lens and the ratio of the trans- 
verse velocity to the minimum impact parameter. Based 
on Monte Carlo simulations of synthetic data sets we give 
the expected accuracies in the parameter determination 
for various masses. 
^ I We argue that pulsars found behind the very center 
1^ of our Galaxy would provide an excellent opportunity to 
test the mass distribution in the Galactic nucleus, and the- 
refore to distinguish between a super massive black hole 
f-^ ^(SMBH) and a super-dense star cluster (SDSC) within the 
I central 0.1 pc. 
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1. Introduction 

Star formation rates increase towards the GC (e.g. Bier- 
mann 1978). Therefore, one can expect a large number 
of massive compact objects, such as neutron stars and, 
perhaps black holes, in central parts of the Galaxy. One 
should also expect a large number of radio pulsars in 
this region, although the present known pulsar popula- 
tion does not show this pattern. In fact, when one corrects 



the known sample of pulsars with respect to various se- 
lection effects, then the largest density of pulsars is found 
near the GC (Lyne et al. 1985, Narayan 1987). There may 
be as many as 10® pulsars in the vicinity of the actual 
GC (Lyne et al. 1985). The most important selection ef- 
fects which bias the sample of discovered pulsars towards 
nearby sources include background radiation, dispersion 
smearing, and scattering pulse broadening. These effects 
can be largely reduced by choosing high observing radio 
frequencies. Although energy spectra of radio pulsars are 
dominated by the low frequency regime, pulsars also seem 
to be quite luminous at high frequencies. Recent detection 
of a number of pulsars at 35 GHz at the Effelsberg Ob- 
servatory indicates the importance of the high frequency 
regime (Wielebinski et al. 1993). It seems that a receiver 
frequency of about 5 GHz should be a reasonable com- 
promise between the lower fiux densities and the largely 
reduced background radiation, dispersion smearing and 
scatter broadening. This allows a large receiver bandwidth 
of about 500 MHz. Detailed considerations for a 100 meter 
radiotelescope equipped with a low-noise receiver (Fig. 1) 
show that such a high frequency search would be sensitive 
to luminous pulsars near the Galactic Center with periods 
longer than about 100 ms. 

Three searches have been undertaken so far at fre- 
quencies near 1.5 GHz. A total of 147 new pulsars have 
been discovered, many of them young and with relatively 
short periods (Johnston 1990, Clifton & Lyne 1986). The 
5 GHz search can avoid scattering broadening for Disper- 
sion Measures (DM) in excess of 1000 cm~^pc. Since the 
base sensitivity for long period signals at 5 GHz is about 
0.2 mjy (Fig. 1), as compared with ~ 1 mjy in previous 
searches at 1.5 GHz (Clifton & Lyne 1986, Johnston 1990), 
one can expect to find new pulsars with relatively fiat spec- 
tra S (X Approximately 50% of pulsars have such a 
fiat spectra with value of a less than 1.7 and the median 
value of a for pulsars detected by Johnston (1990) at 1.5 
GHz is 1.0. Another high frequency search for GC pulsars 
at 2.6 GHz is planned for summer 1995 at Parkes radio- 
observatory (Lyne 1995). 
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Fig. 1. Flux minimum sensitivity curves for a 5 GHz receiver 
with the system temperature Tsys = 30 K installed at a 100 
m radiotelescope. The curves are calculated for DM exceeding 
1000 cm-^pc for both cold (T = 10 K) and hot (T = 100 
K) sky. Total receiver bandwidth assumed is 500 MHz. The 
horizontal line at 1 mJy corresponds to the sensitivity obtained 
in previous high frequency pulsar surveys. 



Small dispersion smearing and negligible scattering at 
5 GHz proves to be crucial in probing deeper into the GC. 
The 5 GHz pulsar search, if conducted with one of the 
world largest radiotelescopes, has the potential of discove- 
ring exotic objects such as newly born pulsars, relativistic 
pulsar neutron-star or pulsar black-hole binary systems. 
The number of such systems in the GC may be large. By 
accurate measurements of the arrival times of the pulses 
one can determine the mass of the pulsar and of its com- 
panion, i.e. if two post-Keplerian parameters are measu- 
red the masses of pulsar and companion are known up to 
the Doppler-factor which is negligible in most cases. This 
may lead to the discovery of the first black hole (see Taylor 
(1993), Taylor & Weisberg (1989) for more information on 
the timing of binary pulsars). Below we discuss the pos- 
sibility of discovering massive objects in the GC even if 
they do not occur in binaries with radio pulsars. 

When the star radiation passes near a massive object, 
its intensity is magnified due to the microlensing phenome- 
non. Several such microlensing effects have recently been 
observed in the direction of the Baade Window in the Ga- 
lactic Bulge (Udalski et al. 1994 1994, Kiraga & Paczynski 
1994). Similar effect should also apply to pulsar signals. 
However, because of the pulsating nature of pulsar radia- 
tion, the fiux magnification will be accompanied by a time 



kova & Doroshenko 1994). In this paper we come to the 
conclusion that regularly timed pulsars in the GC can lead 
to the detection of a number of microlensing events. In the 
case of lenses of more than a solar mass this can lead to a 
very precise mass determination. 

The motion of the stars and gas clouds surrounding 
the very center of our Galaxy is most easily explained by 
a central mass of about IO^Mq within a radius of 0.1 pc 
(Genzel et al. 1994). The observations of a compact, non- 
thermal, radio source termed Sgr A* gives a circumstan- 
tial evidence that our Galactic nucleus houses a SMBH. 
Only recently Narayan et al. (1995) published a model for 
Sgr A* in which a black hole of mass 7 X IO^Mq accretes 
matter with a rate of (1.2 x lO^)aMQ /yi, where a is a 
parameter which describes the viscosity of the gas. The 
model fits the observed spectrum of Sgr A* from radio 
to hard X-ray wavelengths. The observation of a pulsar 
located just behind the Galactic nucleus would provide a 
unique possibility to identify a SMBH and determine its 
mass with high accuracy. Thus the discovery of pulsars in 
the GC could resolve the SMBH-SDSC controversy con- 
cerning the very center of our Galaxy. 

The paper is organized as following: In Sect. 2 we sum- 
marize results concerning the fiux-magnification and an- 
gular separation caused by a Schwarzschild lens. In Sect. 
3 we calculate the time varying travel-time delay caused 
by a passing Schwarzschild lens. We distinguish between 
a weak and a strong travel-time delay event. In Sect. 4 we 
investigate the change of the observed spin parameters of 
the pulsar (period, period derivative, . . .) due to a time- 
varying travel-time delay. In Sect. 5 we present expected 
accuracies of the parameter determination in the case of 
weak travel-time delay. In Sect. 6 we present expected ac- 
curacies in the parameter determination in case of strong 
travel-time delay. In Sect. 7 we compare the travel-time 
delay caused by a SMBH and the travel-time delay caused 
by a SDSC in the Galactic nucleus. In Sect. 8 we present 
statistical investigations. We calculate the probability of 
observing a time-delay effect toward the GC and the pos- 
sibility of finding a pulsar behind a possible SMBH in the 
Galactic nucleus. Finally, we present our conclusions in 
Sect. 9. 

2. Flux magnification and angular separation 

Let us consider the infiuence of a compact astronomical 
object with a mass M on the radiation of a pulsar, when 
it passes close to the pulsar-Earth axis. This could be an 
ordinary star, a white dwarf, a neutron star or a black hole 
and, generally, we shall call it the lens. The corresponding 
lensing geometry for such a Schwarzschild lens is presented 
in Fig. 2. PSR denotes the position of the pulsar and SSB 
denotes the position of the barycentre of the solar system, 
by what we mean the proper reference frame attached to 
the mass centre of the Solar System. 
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Fig. 2. Lensing geometry of the PSR-M-SSB system. Here 
h is the impact parameter and {-|- are the distances in the 
lens plain between the mass and the two different light rays 
connecting PSR and SSB. 



As for ordinary stars the pulsar radiation will experi- 
ence a flux magnification. If the observed intensity in case 
of no lens is equal to 7o then the observed intensities of 
the two images, 7_|_ and 7_ , obey the following relations: 
(Refsdal 1964) 
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where / is a dimensionless impact parameter 
/ = h/RE 
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is the Emsiem radius. Here Rs is the Schwarzs child radius 
of the lens 
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The ratio of the observed intensities of the individual 
images is 
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If the impact parameter h is much larger than the Ein- 
stein radius Re, then the contribution of the second image 
to the total brightness is completely negligible. For in- 
stance if / > 2.85 then the contribution of the second 
image is less than 1%. 

In the case when both pulsar and lens are located in 



the GC region we find Dmp ^ Dm 
Einstein radius is approximated by 



D„ . Therefore the 
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For normal stars acting as lenses in the GC the Einstein 
radius is a few AU. For a SMBH of ~ 10^ Mq we find an 
Einstein radius of the order of 10^ AU. 

The angular separation of the two images is (see Refs- 
dal 1964) 



Ae = ^VPT-i. 
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If pulsar and lens are located in the GC we have Dm ~ 8 
kpc. Using Eq. (7) we find 



M « 3'.'6 X 10-4/f^ ( — 



1 kpc 



(9) 



For a SMBH (M ~ IO^Mq) we expect M to be of the 
order of arc seconds. 

3. travel-time delay 

The pulsating nature of pulsar radiation offers the possi- 
bility of observing another aspect of gravitational lensing, 
namely the change in the light travel-time. There are two 
effects that contribute to the light-travel time of the pul- 
sar signals (see Cooke & Kantowski 1975, Schneider et al. 
1993). First, the defiection of the light ray increases the 
length of the light path {geometrical time delay). Secondly, 
the light traverses the gravitational field of the lens and 
therefore suffers a potential time delay. If & ^ Dmp, Dm, 
the light-travel time for a Schwarzschild lens is given by 
(see Appendix A) 



T± = TM {4:F+ ^ - 2 In T'i) -I- const. 



where we introduced the "Schwarzschild time" 



TM = Rs/c. 



(10) 
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Since the lens is moving, the light-travel time depends 
on the time of observation in a characteristic way. Writing 
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where the definition 



we find that 
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/ = /m Vl + s2 , s = q(t-To), 
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Here To is the time of the closest approach of the lens to 
the PSR-SSB line. 

In the case of a weak lensmg event (fm ^ 1) we can 
neglect the geometrical time delay in Eq.(lO) for the first 
light ray (+). The time delay is then dominated by the 
potential time delay 

r_|_ = —tm ln(l + s^) + const. (15) 

In the solar system this is known as the Shapiro delay 
(Shapiro 1964) which is tested with a measurement preci- 
sion of 0.1% (see e.g. Will 1992). According to Eq. (1) the 
intensity of the second image (— ) is negligible. 

In order to measure the travel-time delays, we need a 
number of precise measurements of topocentric pulse arri- 
val times (TOAs). Let the parameter s = sq corresponds 
to the beginning of TOA monitoring at the initial time 
t = to- We define 



Ar = r(/)-r(/o). 



(16) 



For the brighter image (-|-) the maximum delay occurs 
at s = 0, corresponding to the closest approach of the lens 
to the line-of-sight when h = hm 
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The time between the beginning of regular observati- 
ons {to) and the maximum (To) written in more suitable 
units is 

hm[AU] 



To-to = -(4.7 yr) sq 



-(4.7yr)so Q"'- (18) 



v±_ [km/s] 

Let us consider lensing objects with stellar masses 
M = IMq . ..20Mq. For ^ 100 AU, v_l ^ 300 km/s, 
and So = —1 we find that 



To - to ~ 1-6 yr. 



(19) 



and the corresponding maximum differential Shapiro de- 
lay is 
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A more detailed investigation of Eq. (16) for the case 
of a weak travel-time delay (cf. Eq. (15)) is presented in 
Fig. 3. 

The actual GC (Sagittarius A West) is likely to con- 
tain a SMBH with a mass M ~ IO^Mq. Radio-continuum 
VLBI observations show a core source not larger than se- 
veral AU in size. One can therefore expect that radiation 
of the innermost GC pulsars will be gravitationally af- 
fected by the putative SMBH. The Shapiro delay will be 
strongest for pulsars with small angular separation from 
the Sag A West. The fact that in this case the lens is at 
rest and the pulsar is moving does not change the con- 
siderations given above. As can be seen from Fig. 3 the 
Shapiro delay for M/Mq ~ 10® is of the order of several 
seconds. 




Fig. 3. Differential Shapiro delay [in 1O"''(M/M0) s] as a func- 
tion of time t — To [in years]. 



4. travel-time delay induced period variations 

In this section we investigate the infiuence of the travel- 
time delay on the observed period P and the observed 
period change P of the pulsar. Explicit equations and fi- 
gures we shall give only for the weak travel-time delay, 
see Eq. (15). The infiuence of a strong travel-time delay is 
qualitatively the same. 

Let us denote the observed period and the observed 
period change in the case of no lens ("intrinsic") by Pi 
and Pi, respectively. To simplify the notation we shall use 
the definition 



M = M/Mq . 



(21) 



The difference in arrival time of two consecutive pulses 
is the observed period of the pulsar, which is given by 

dr 

P = t2-ti = Pi + T{t2) - T(ti) ~ Pi + Pi-^ih). (22) 



Since ds/dt = q (Eq. 14) we find 
2s 



T+ ~ Ti 1 - TMq 
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One can see that P is equal to Pi for s = and s = ±cxd 
and the deviation of P from Pi is maximum for s = ±1. 
Thus 



Pi 
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TMq ~ 6.58 X lO-^^'A^Q. 



(24) 



For a weak lensing event the change of the period _P_|_ with 
time is shown in Fig. 4. 




The time derivative of the observed period is 
,2(1 -.2) 



P+ ~ Pi - PiTMq' 
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The change of the period due to the travel-time delay is 
maximum for s = 0: 
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(8.8 X 10-22 s-1)A^Q2 



For example, if t;j_ = 300 km/s, h„ 
1 ... 20 Mq we find that 
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Variations of _P_|_ with time for a weak lensing event are 
presented in Fig. 5. 

The measurement precision for P / P for a typical mil- 
lisecond pulsar observed for several years is in the range 
of 10"^^ . . . 10~^® s~^. Therefore a lensing event in our 
Galaxy with ~ 100 AU might just be observable in the 
case of a stellar mass. 

The situation is completely different for a central 
SMBH with M ~ IO^Mq, since the mass is scaled by 
factor of 10®. In this case the modulation of pulsar period 
by microlensing is comparable with a slow-down rate due 
to the magnetic dipole radiation, which is typically 10~^® 



Finally, we would like to mention that the time varying 
travel-time delay can infiuence significantly the observed 
second period derivative _P_|_. If ^ 1 than _P_|_ is maxi- 
mum for s = ±(-^2 — 1) and 
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(8.6 X 10"^° s-2)A^q3 



(28) 



(26) Thus, even a rather weak travel-time delay (&„ 
100 AU,!"! ~ 100 km/s) could cause a "braking index" 
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that is much bigger than 3 (ris = 3 is expected for rotatio- 
nal energy loss by dipole radiation). For example, in case 
of PSR B1937-I-21 the observed braking index is about 
4500. This can be explained by a travel-time delay caused 
by a lens of a few solar masses (see Wex 1995). A number 
of pulsars show the value of considerably larger than 
3. Perhaps this is a result of the travel-time delay caused 
by a massive object passing near the line-of-sight. 

5. Parameter analysis for a weak lensing event 

For lenses with stellar masses {M = IMq ■ ■ ■20Mq) the 
Einstein radius is expected to be of the order of a few AU, 
as shown in Sect. 2. Thus in this case a strong lensing of 
the pulsar radiation (/ <, 1) is very unlikely (cf. Sect. 8). 
Therefore, we will now concentrate solely on weak lensing 



events to the next section. 

In this section we investigate the expected accuracies 
for the parameter determination in the case of a weak 
lensing event (/„ ^ 1). Since the second image (— ) is 
much fainter than the first one (+) there is only one pulsar 
signal for timing, which practically does not change in 
intensity, i.e. the change in intensity is too weak to allow 
for some parameter determination. Since the arrival time 
of pulsar signals can be measured with very high accuracy 
we have a chance to obtain information about the lens 
by observing the travel-time delay caused by the lensing 
mass. We present a ttmmg- formula for pulsars that show 
a weak travel-time delay in their TO As. 

Let us denote the epoch of observation by . A weak 
travel-time delay is characterized by three parameters, M , 
q, and To (see Eq. 15). We denote the "intrinsic" frequency 
and its derivatives with respect to time corresponding to 
the barycentric arrival time t = t((, hy i>, i>, i), . . . (Note: 
"intrinsic" denotes the frequency and its derivatives mea- 
sured in the absence of any lensing mass.) 

The "intrinsic" frequency corresponding to a certain 
barycentric arrival time t is 



l^i = ly + v{t - t^)) + \v{t - tQ)f + . . . 

The observed frequency at the time t is 
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Since dr/dt ~ 10 ^^AiQ we find (neglecting i) and higher 
derivatives of i>) 
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The pulsar phase at the barycentric arrival time t is given 

by 



(33) 



i>(t)rft = $(t;^Ao,J^>,M,g,To). 
Thus 

$(t;<?io,J^>,M,g,To) = 

<?io + v{i - i®) + \v{i - i®? + J^Ash, 
where 

Ash = -rMln[l + g2(t-To)2] . 



(34) 



(35) 



The inverse function of $ ($ ^) gives the barycentric 
time at which the pulsar shows the phase : 



t = $-i(<?i;<?io,J^>,i>,M,g,To). 



(36) 



Having N measurements {4>k,tk), where 4>k and are 
pulsar phase and the corresponding (barycentric) arrival 



6i . . .6„ (n = number of parameters), by the least square 
fit to the data: 



N 



k = l 



Pi 



Min. 



(37) 



where is the error in the barycentric arrival time of the 
k^^ measurement. It is convenient to use $ instead of $~^: 
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We end this section with an investigation of the expec- 
ted accuracies for derived parameters. We will introduce 
two length scales denoted by TZi and TZ'z- Here TZi is the 
maximum value of the impact parameter that allows a 
reliable parameter determination and by TZ'z we denote a 
value of the impact parameter at which the travel-time 
delay is just visible in the residuals of a (/)o-J^-i'-fit. 

For the numerical investigations we assume a 10 years 
observation of a 100 Hz pulsar where TOAs are taken twice 
a month. Moreover we assume t;j_ = 300 km/s. To — to = 5 
yr, i> = 0, and Re (see Eq. 3). W^e further assume 

that the pulse profile is sampled with 1024 bin and that 
the l(T-error for timing is 0.1 bin. 

Monte Carlo simulations of synthetic data sets lead to 
the values given in Tables 1. For lenses below 0.5 Mq the 



Table 1. Weak lensing event 



M 


0.1 1 10 20 


bm [AU] 
AM/M 

Aq/q 

A(tffi-To) 
(tffi-To) 


— 10 50 100 

— 10% 5% 7% 

— 20% 5% 5% 

— 2% 0.5% 0.5% 


7^l [AU] 
7^2 [AU] 


— 30 100 200 
20 150 350 500 



parameter determination is rather poor. But there is still 
the chance to see the effect in the residuals. We find as a 
good approximation 



r250A^i 1 if 0.1 < < 0.5, 
7^2[AU] « <^ 

[l50A^°^ if 0.5 < < 20. 



(39) 



This result for TZ2 was obtained for v±_ = 300 km/s, 
though for M <, 20Mq TZ'z proves to be not very sensi- 
tive to Vi_ as long as t;j_ ^ 100. 

If the impact parameter exceeds TZ'z then the travel- 
time delay will not show up in the residuals. Generally 
speaking, this is a situation where q is very small, i.e. 
is very big and/or t;j_ is very small. In this case Eq. (35) 
reduces to 
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certain time t is given by 

^0 = (j>o - VTMq^(t(s -Tof, 
V = v[l-2TMqH-t®-To)\ , 
v = v — 2vtm({'^ ■ 



(41) 



Therefore the travel-time delay is hidden in the observed 
period and period derivative, i.e. it is not possible to ex- 
tract any information on M , q, and To from the TO As wi- 
thout making assumptions about the intrinsic frequency 
(i>) and the intrinsic frequency change (i>) of the pulsar. 

6. Parameter analysis for a strong lensing event 

(fm <> l) 

If the impact parameter h is of the order of the Einstein 
radius Re we call this a strong microlensmg event. In this 
case the second image is of comparable brightness, see Eq. 
(5), e.g. for / = 0.5 the second image (— ) has approxima- 
tely 40% of the intensity of the first image (-|-). Receiving 
two pulsar signals of the same pulsar can be used to ex- 
tract more efficiently information from the observations. 

The radiotelescope will receive two pulsar signals with 
different periods and different intensities. For the phase 
between the strong pulse (-|-) and the weak pulse (— ) we 
find 



6 = — mod 1, 
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r_ - r+ = TM 



/V/2+4 + 21n 



(42) 



(43) 



see Eq. (10). (Eq. (43) was obtained in a more general 
form by Cooke & Kantowski 1975). Using the relation 



T+ = Tm 



r — 1 



In 'i 



(44) 



(see Krauss & Small 1991) we get (5 as a function of r. In 
principle this could be used to determine the mass M if 

P+>T.-T+. 

6.1. Stellar masses 

For stellar masses {M = 1 . . . 2OM0) tm is of the order 10 
jjLS to 0.2 ms. But the minimum 10%-width of the pulse 
profile at 4.75 GHZ is (see Kramer et al. 1994) 



1^10% « 9° X {P[s])-° \ 

which corresponds to a time width of 

,5tio%«(25ms)x(T[s])°5. 



(45) 



(46) 



(5^10% and so one will not be able to see two separate pul- 
ses, and the measurement of r_ — r_|_ will not be possible. 
There is only one pulse brightened by the factor A given by 
Eq. (1). Only in the case of a millisecond pulsar {v ^ 100 
Hz) and M ^ 2OM0 we could expect a two-pulse profile. 

As in the case of microlensed ordinary stars the change 
of A with time can be used to extract some information 
from the observations. The problem here is that pulsar 
radiations shows strong luminosity variations caused by 
scintillation (see Ricket 1977), intrinsic instabilities, etc., 
which can limit the accuracy of parameter determination 
quite strongly. 

As in Sect. 5 the most promising method to get infor- 
mation on the lens is the travel-time delay. But this time 
we have to use Eq. (10) instead of the simple Eq. (15). We 
get fm as a further parameter to characterize the strong 
microlensing event. Analogously to Sect. 5 we find 



^{t;(f>o,i^,0,M,frn,q,To) = 
where 

\+ = tm{^F^^ -'2\nF+) , 



(47) 



(48) 



Similarly to Sect. 5 we performed Monte Carlo si- 
mulations of synthetic data sets to estimate the expec- 
ted accuracies in the parameter determination. We arri- 
ved at the conclusion that for stellar masses it will not 
be possible to get a reliable fit for all seven parameters 
((/)o, I', v, M, fm, q, To) unless timing accuracy is increased 
by at least an order of magnitude. Especially the para- 
meters fm and q appear to have rather large errors. The 
most practicable method for M ~ IOMq and fm ^ 0.5 
is to limit fm by the determination of the maximum fiux 
magnification A^nax = A(fm), see Eqs. (1), and then fit 
for the other parameters. Table 2 shows the result for 
a 5 years weekly observation of a 100 Hz pulsar around 
the maximum of the microlensing event (v = — 10~^®s~^, 
v±_ = 300 km/s). We assumed an accuracy of 20% for the 
measurement of the maximum fiux magnification A^nax- 



Table 2. Strong microlensing event: M = WMq, Re = 10 AU 





6„ = 2 AU 6„ = 5 AU 


fm 


0.17. ..0.25 0.40... 0.66 


AM/M 

Aq/q 

A(tffi-To) 
(tffi-To) 


1% 2% 
16% 13% 
0.2% 0.2% 



In case of a sun-like mass the strong travel-time delay 
will clearly show up in the residuals. But even if we use 
Amax to extract fm from the observations, there will be no 



we can restrict the error in to about 20%, Monte Carlo 
simulations show that there will be nearly no restrictions 
on the parameter q. If we restrict q by some plausible 
a priori assumptions, e.g. about v_i, we can still extract 
acceptable values for M and To from the data with an 
error that does depend only very little on the error made 
in the a priori assumption. 

6.2. A SMBH m the Galactic Center 

For a SMBH of IO^Mq hosted in the very center of the 
Galaxy we find the Einstein radius to be several orders 
bigger than for stellar masses 



Re ^ (2900 AU) 



Dmp 
1 kpc 



1/2 



(49) 



But a pulsar located a few 100 pc behind the GC and ta- 
king part in the Galactic rotation will have a transverse 
velocity of v±_ ~ 150 km/s (see Kent 1992). Therefore the 
pulsar will cross in the lens plane a path of 32 AU in one 
year, which is by far smaller than the Einstein radius. The- 
refore / will change only by a small amount A/ during an 
observation of a few years and we need / ^ 1 to see a no- 
ticeable change during such an observation. Thus we can 
expect to see both images where 1 < r ^ 10. Recording 
a typical "bell-shaped" curve (see Fig. 3) in the timing 
data would require an observation of about a 100 years. 
In Table 3 we present the changes of various fiux and time 
parameters during the last 2 years before reaching the mi- 
nimum of the impact parameter (t = To — 2 yr . . . To), for 
a pulsar located a 100 pc behind the SMBH. 



some kind of distinct outbursts in the pulse emission. 

The first two consequences imply that there is only 
very little chance to extract information from the measu- 
red intensities of the pulsar signals. For ^ 1000 AU we 
can use Eqs. (5) to extract a crude estimation for / from 
the observed intensities. The corresponding error is 



A/ 



VP 



4 Ar 

r 



(50) 



Much more promising is the use of Eq. (43). Here we 
have to measure the phase 6 as a function of time. Figure 6 
illustrates the change of the observed pulse profile within 
a period of two years before To. 
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Table 3. SMBH: M = 


IO'^Mq, 


Re = 


1000 AU 




[AU] 


/*+ 






r 




500 


1.58. . . 1.59 


0.58. . 


0.59 


2.71. . 


2.69 


200 


2. 93... 3. 04 


1.93. . 


2.04 


1.52. . 


1.49 


50 


6.72. . . 10.5 


5.72. . 


9.52 


1.17. . 


1.11 


6™ [AU] 


Ar+ 


Ar 




T- — 




500 


61 ms 


-101 


ms 


20.04. . . 


19.88 s 


200 


173 ms 


-213 


ms 


8.27. . . 


7.88 s 


50 


583 ms 


-623 


ms 


3.18. . . 


1.97 s 



The following three consequences can be extracted from 
Table 3: 

— There is practically no variation in the ratio of the in- 
tensities of the two images within years of observation 
for all values of We can safely use r ~ const. 

— Only for small values of / there is a measurable change 
in //_|_ and . 

— r_ — r_|_ is in general much greater than a typical period 
of a pulsar. Therefore it is not possible to read off 
r_ — T_|_ from the pulse profile. Consequently we cannot 



phase 



Fig. 6. Change of the pulse profile due to gravitational lensing 
by a SMBH (M = IO'^Mq) in the GC. The pulsar has a period 
of 0.25 s and is located 100 pc behind the GC. 6„ = 500 AU, 
v± = 150 km/s. The pulse profiles are shown in steps of 0.5 
years where the last pulse profile corresponds to the time To. 



Monte Carlo simulations of synthetic data sets show 
that by fitting 



P+6(t) = T+-T_-NsP+, 



mem 



(51) 



to the observations we have to distinguish three cases: 

I) Eq. (43) can be approximated by 

P+6(t) = 2TMfmq^(t - Tof + (4:TMfm " mP+). (52) 

In general, has to be a few times larger than Re. 
But then ^ 1 and therefore in practice we will 
not be able to see the second image. Eq. (52) does not 
allow mass determination. 



to fit all five parameters (M, fm, (l,To, Ns). Here we 
use Eq. (5) to get a crude estimation for the value 
of fm = f(To)- (The time To can be extracted quite 
easily from the observation data with sufficient high 
accuracy.) We will use this result in Eq. (51) and then 
fit for the parameters M, q,To, Ns . An error of 25% 
in fm leads typically to an error of a 10. . .20% in the 
determination of the mass M . 
Ill) h <, 500 AU. We have to use the full Eq. (43) to fit 
the observed data. Table 4 shows the expected accu- 
racy for the mass determination for a 10-years monthly 
observation of a 250 ms pulsar, to — Tq = —5 yr. 
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Table 4. SMBH: M = 10^ Mq, Re = 1000 AU 
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100 


50 
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Fig. 7. Residuals of a i/iq — i' — I'-fit for a 10 years observation 
of a pulsar that is lensed by a IO'^Mq SMBH in the GC where 
bm = 2000 AU, vi_ = 150 km/s. Re = 1000 AU, and To -to = 5 

yr- 



Concerning cases II and III it is important to note that on 
one hand the period of the pulsar is very well known, and 
on the other hand the quantity Ns is an integer number. 
In case of long-period pulsar Ns can be determined with 
rather high confidence, and thus it can be advantageous to 
use Eqs. (51) and (44) to determine the mass, especially 
in case II. 

As in the case of stellar masses we can also use the 
TOAs of the brighter pulse (-|-) to do a parameter analysis 
by fitting Eq. (47). Monte Carlo simulations of synthetic 
data sets show that this method is practicable for ^ 
300 AU . The expected accuracies are worse by more than 
a factor of two compared with the accuracies achieved in 
the previous method, i.e. by fitting 6(t). 

The infiuence of the travel-time delay can be seen in 
the residuals of a 10 years observation of a 10 Hz pulsar 
up to an impact parameter of about 2000 AU (see Fig. 
7), and thus 7^2 ^ 2000 AU for a IO^Mq SMBH in the 
Galactic center. Only for ^ 2000 AU fitting Eq. (47) 
with tm = would provide a good fit. 

7. super-massive black hole versus super dense 
star cluster 

So far we have payed attention only to the Schwarzschild 
masses as candidates for lenses. Since at the present time 
it is not clear whether the very center of our Galaxy hosts 
a SMBH or a super-dense star cluster (SDSC) we will 
investigate the possibility of distinguishing between these 
two models by observing a lensed pulsar in direction of 
the Galactic nucleus. 

As mentioned in the Introduction, observations in- 
dicate that a mass of IO^Mq is resident within a 0.1 
pc (~ 20 000 AU) radius around the very center of the 
Galaxy. This corresponds to a mass density of about 



IO^Mq/pc^. We use Plummer's law (see Binney & Tre- 
maine 1987) 



g(R) 



[1 + (i?/i?*)2]5/2 



(53) 



to model the density distribution of the putative SDSC. 
Here go is the central mass density and i?* is a length 
that characterizes the size of the SDSC. Let us denote the 
total mass of the star cluster by . Integration of Eq. 
(53) leads to 



M, = ^TTpoRt 



(54) 



About 85% of Mh, are within a radius of 3/?*. The corre- 
sponding surface mass density is 



where 



So 



[1 + (i/R*) 



212 ' 



So - |po-R* - 



i Ml 



(55) 



(56) 



is the central surface mass density. If So is smaller than 
the critical surface mass density 



AttG DmD 



Mp 



(57) 



we will observe only one image of the pulsar. If So > S^r 
we will have one or three images, depending on h (see 
Schneider et al. 1991). For a lensing event in the GC (Dp 
Dm) we find 



So/Scr ^ 1.9 X 10" 



10^'MqJ vipc 



D 



Mp 



1 pc 



(58) 



well below 1 for Dmp ^ 1 kpc and thus we expect only 
one image for all impact parameters h. 

To shorten notation and facilitate the comparison bet- 
ween the SMBH and the SDSC we define 



Rs 



Re = \ / 2i?s 



Dm Dmp 



We define the dimensionless parameter 

Ko = So/Scr = Rl/Rl (60) 

For small impact parameters (h R*) the lens equa- 
tion for Plummer's model leads to (see Appendix A) 



1 - Ko 

The corresponding travel-time delay is (see Appendix A) 
T = —tm, (f^o^ — ^) J + const., (62) 
where we used the definitions 

TM,=Rs/c, f = b/RE. (63) 
Since 

/ = L Vl + s\ .s = q(t-To) (64) 
we obtain 

T = —C^,(t — Tq)^ + const., 



TM, / vj_ 



(65) 



1 - Ko \Re, 

Thus, in contrast to the SMBH, 
would fit the observed TOAs. We find 

<j> = ^o + v(t- t(s) + ^v(t - ^e)^ 

^0 = (j>o- vC^{t(^ - Tof, 
V = v[l-2CS®-To)], 
v = i/ — 2i>C^,. 



(66) 



(67) 



For a SDSC with = IO^Mq, i?* = 0.1 pc. Re = 1000 
AU, and v±_ = 150 km/s one finds 



i>-^ = iy(2.0x 10-^^s-i). 



(68) 



Since v ~ —10~^^ for most pulsars, one does not expect 
to discover the infiuence of a SDSC in the TOAs. Note: 
The infiuence of the SDSC is independent from . 



Based on the numerical investigations in Sect. 5 and Sect. 
6 we give some coarse estimations for the probability of mi- 
crolensing of pulsar radiation in the GC. (See also Krauss 
& Small 1991). 



(59) 

8.1. Luminous matter 



To estimate the lensing probability by ordinary stars we 
shall calculate the expected number N of stars between 
the pulsar and the Earth which are able to infiuence no- 
ticeably the TOAs. 

For the density of luminous sources in the Galactic 
plane we use an empirical model of Bahcall & Soniera 
(1980): 



(61) Qir) 



Qdisk{r) + Qbu\ee{r) 

7 Mq 



A X 10' 



kpc-- 



exp 







3 X 10^ 



M, 







kpc-- 



1 kpc 



3.5 kpc 

-1.8 

exp 



(69) 



1 kpc 



Here gdisk and £>buige are the contributions of the disk and 
the bulge, respectively, Rq is the distance of the Sun from 
the Galactic Center (Ri 8 kpc) and r is the distance from 
the Galactic Center. £'buige('') is not valid for the very 
central region (r ^ 1 pc). 

The number density of the stars that have a mass bet- 
ween A4 and A4 + dAi can be estimated using the initial 
mass function given in Kroupa (1995). We find 



£,{M,r)dM 



if 0.08 < < 0.5, 



1.4 



(m^^) ''■^ 1 0-5' if 0.5 < < 1.0, (70) 



0.52 2A^-2 ^ if 1.0 < < CXD. 



Thus the number of stars that infiuence noticeably the 
TOAs is 



N 



T:[n2{M)fS,{M,l) dM 



(71) 



0.08 



where / parameterizes the pulsar-Earth distance. The den- 
sity Q is given by Eq. (69) and the radius TZ'z we take from 
Sect. 5, Eq. (39). Integration in Eq. (71) over the mass Ai 
leads to 



# = 1.6 X 10- 



/.£)Q[kpc] 

-12/ e(O[M0/kpc-']rf/. 
■Jo 



(72) 



Figure 8 shows the resulting number N for a pulsar 
located at r = 0.5 kpc as function of the angle between 
the directions GC-SSB and GC-PSR. Strictly speaking. 
Fig. 8 is valid for the observation of a 100 Hz pulsar where 




180 



Fig. 8. Expected number of lenses N for a 100 Hz pulsar 
at 0.5 kpc distance from the GC as function of the angle a 
between the directions GC-SSB and GC-PSR. v± ~ 300 km/s 
is assumed. 



v±_ ~ 300 km/s was assumed for all the lenses. For the 
majority of the pulsars, which has frequencies well below 
100 Hz, the number N is smaller by an order of magnitude. 

Figure 8 suggests that for pulsars located behind the 
GC there is a high probability that the timing data are 
influenced by some travel-time delay. If a considerable part 
of the dark matter in our Galaxy exists in form of massive 
objects (M ^ O.5M0) the lensing probability is increased 
by an order of magnitude. 



and, perhaps, even a SMBH. In this paper we consider the 
relativistic interaction of pulsar radiation with a moving 
massive object. Based on extensive Monte Carlo simula- 
tions of synthetic data sets we argue, that the important 
new effect to be taken into account in pulsar timing is the 
travel-time delay. This could lead to a very good deter- 
mination of the mass of the lens as it was shown in Sect. 
5 and Sect. 6. The statistical estimates of Sect. 8 imply 
that the probability for observing a travel-time event with 
a lens mass between 1 and 20 Mq is quite reasonable due 
to the high density of stars in the galactic bulge and the 
large radius of influence TZ2 even for small masses. 

In Sect. 6 it was shown that the travel-time delay could 
provide a unique possibility for the detection of a SMBH in 
the GC. Regular timing observations over a period of a few 
years would lead to a mass determination with an accuracy 
of a few percent. As calculated in Sect. 7 the travel-time 
delay caused by a SDSC in the Galactic nucleus will lead 
to TOAs which are indistinguishable from the TOAs of an 
unlensed isolated pulsar and thus they are quite different 
from those of a SMBH. Thus the detection of a lensed 
pulsar in the very center of the Galaxy could bring to 
the end the SMBH-SDSC-controversy. Unfortunately the 
situation having a pulsar within a radius of a 1000 AU 
around the galactic nucleus (Sag A*) seems to be not very 
likely. 
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8.2. The SMBH model 

Finally we would like to estimate the expected number of 
pulsars behind a putative SMBH in the GC. For a SMBH 
of IO^Mq we have 7^l ~ 1000 AU (7^2 ~ 2000 AU). Using 
Eq. (69) we find a mass of 70 Mq for 7^l (280 Mq for TZ-j) 
between 1 pc and 1 kpc behind the SMBH. If we assume 
that 0.1% of the mass is provided by neutron stars, then 
there should be approximately 0.07 (0.3) neutron stars be- 
hind the SMBH. If 10% of these neutron stars are visible as 
pulsars, we can expect more than 0.007 (0.03) pulsars suf- 
fering an observable Shapiro delay caused by the SMBH. 
This number is further decreased by the fact that one ex- 
pects to find mostly the young pulsars in the GC, which 
are both bright and have wider beams. 

9. Conclusions 



A. The travel-time delay 

By travel-time delay of a light ray we denote the (coordi- 
nate) time delay of this light ray relative to the undefiected 
light ray propagating in an Euclidean background metric. 
The travel-time delay for the light ray that intersects the 
lens plane at the position ^ is given by (see e.g. Schneider 
et al. 1991) 



cr = ^^^[mf-m + const. 



(Al) 



where the defiection potential tp and the defiection angle 
a are given by 



^(0 = ^ /rfYS(nin^'^ 



(A2) 



The discovery of a large number of pulsars in the GC seems 
to be only a matter of time. The GC region contains pro- 



«(0 = V^(0, (A3) 



density of the deflector onto the lens plane and is a 
typical length scale in the lens plane. 

For axially symmetric lenses, i.e. S(^) = S(|^|), we 
flnd (see e.g. Schneider et al. 1991) 



const. 



(A4) 



and 

m = ^ ['deem')- 



A.l. The Schwarzs child lens 



If the deflected light ray propagates in a Schwarzschild 
spacetime characterized by the mass M then 



We obtain 



a(0 



Since 



4GM, / ^ 



In 



Re 



AGM 



^ = ^± = ^F±, F± = VP + ^±f (A9) 



(A6) 



(A7) 



(A8) 



we obtain further 

CT± = Rs (4-F±^ - 21n_F±+) + const. 
A. 2. Plummer's model 



(AlO) 



If our deflecting mass is a spherically symmetric star clu- 
ster with a density distribution that follows Plummer's 
law then the corresponding surface mass density is given 

by 



1 Mh. 



(All) 



[l + (e/i?*)2]2' TT 

Using Eqs. (A4),(A5) we obtain 

^.(0 = Rs ln[l + {ilR*f] + const. (A12) 
and 



"(0 = 2i?..^,^^, 
The lens equation is 
Dm Dmp 



(A13) 



(A14) 



i.e. Re < R*, then Eq. (A14) has only one solution. For 
small impact parameters h (^ i?*) we flnd 



Ko 



R% 



(A15) 



1 — Ko ' " R* 

Consequently, using Eq. (Al), we flnd 

CT = -Rs{l- Ko)-^P + const., f = b/RE. (A16) 
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